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ABSTRACT 


This report discusses the system developed to process digitized 
telemetry data from the intensity monitoring spectrometer flown on 
the Orbiting Geophysical Observatory-F (OGO-F) satellite. Functional 
descriptions and operating instructions are included for each program 
in the system. The Automated Production Logging System developed 
and implemented during this report is also discussed. 
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FOREWORD 


The design and implementation of the data processing systems 
herein are the result of analytical research performed for the Solar 
Ultraviolet Branch of the Aeronomy Laboratory, Air Force Cambridge 
Research Laboratories, Air Force Systems Command, Hanscom Field, 
Bedford, Massachusetts. 


PRECEDING PAGE BLANK NOT FILMED 


\ 


ili 



Table of Contents 


Section Page 

Abstract i 

Foreword iii 

List of Illustrations vii 

1 Introduction 1-1 

2 Scope of the Data Processing System 2-1 

3 Experiment Technical Description 3-1 

4 Format of Digitized Data Tapes 4-1 


5 Description of Computer Programs and 

Operating Instructions 

5.1 Programs RDLP and PKGROUP 

5.1.1 Functional Description 

5.1.2 Operating Instructions 

5. 1.2.1 System Control Cards 

5. 1.2. 2 RDLP Fortran Card Deck 

5. 1.2. 3 RDLP Data Cards 

5. 1.2. 4 PKGROUP Fortran Card Deck 

5. 1.2. 5 PKGROUP Data Cards 

5. 1.2. 6 Output from Programs 
RDLP and PKGROUP 

5.2 Background Application Program 

5.3 Peak Intensity List Program 

5.4 Relative Intensity Analysis Program 
(OGOREL) 

5.5 Peak Plot Program 

5.6 Statistical Summary Program (SATAV) 

5.7 Statistical Summary Plot Program (SUMPLOT) 


6 OGO-F Production Master File Logging System 6-1 

6.1 Sources of Data 6-1 

6.1.1 OGO-F Data Tape Packing Sheet 6-1 

6.1.2 OGO-F Data Parameters from 

Research Personnel 6-1 

6.2 Master File Data Record 6-2 


5-1 

5-1 

5-1 

5-8 

5-8 

5-9 

5-9 

5-14 

5-14 

5-15 

5-15 

5-16 

5-18 

5-19 

5-22 

5-24 


TOBKEDiNe PAGE 


Not filmed 


V 



(Cont *d) 


Table of Contents 


Section 


Page 


6,3 OGO-F Master File Printout 

6-2 


6.4 Updating Production Master File 

6.5 OGO-F Chronological Printout - 

6-3 


15-Minute Intervals 

6-3 

7 

OGO-F Interactive Graphics Application 
Program ( OGOGRAF ) 

7-1 

8 

References 

8-1 

9 

Related Contracts and Publications 

9-1 

10 

Contributors 

10-1 


vi 




List of Illustrations 


Figure 


Page No 

4.1 

Data Format for Experimenter 9 

4-4 

5.1 

File Identification Data 

5-27 

5.2 

Long Scan Plot Frame 

5-28 

5.3 

Short Scan Plot Frame 

5-29 

5.4 

Wavelength and Relative Intensity Values 

5-30 

5.5 

RDLP and PKGROUP Deck Setup 

5-32 

5.6 

CDC 6600 Input Card and Plot Request Card 

5-33 

5.7 

Peak Plot Frame 

5-34 

5.8 

Statistical Summary Program Printout (SATAV) 

5-35 

5.9 

Statistical Summary Plot Frame (SUMPLOT) 

5-36 

6.1 

OGO-F Data Tape Packing Sheet 

6-4 

6.2 

Format of OGO-F Master File Data Record 

6-5 

6.3 

Format of Punched Card Generated from Master 
File Data Record 

6-7 

6.4 

OGO-F Master File Printout 

6-8 

6.4.1 

OGO-F Master File Printout 

6-9 

6.5 

Format and Description of Master File Printout 

6-10 

6.6 

Punched Card Format to Update Production 
Master File 

6-11 

6.7 

File Status - Numerical Code Definitions 

6-12 

6.8 

OGO-F Update Printout 

6-13 

6.8.1 

OGO-F Update Printout 

6-14 


vii 



Figure Page No. 


6.9 

Format and Description of Production 
Log Printout 

6-15 

6.10 

OGO-F Chronological Report - 15-minute 
Intervals 

6-16 

7.1 

Sensor 

3600 

Corrected Sensor Data 

7-7 

7.2 

Sensor 

3600 

Uncorrected Sensor Data 

7-8 

7.3 

Sensor 

3600 

Y-Scale Modification 

7-9 

7.4 

Sensor 

3600 

Scissor Application 

7-10 

7.5 

Sensor 

1200 

Before Background Subtraction 

7-11 


7.6 


Sensor 3600 


Redefined Backlash Correction 


7-12 



1.0 INTRODUCTION 


This report describes the data processing system designed to 
analyze digitized data acquired from a spectrometer experiment 
flown on the Orbiting Geophysical Observatory (OGO-F) satellite. 

The report focuses upon the logical procedures and techniques 
developed to perform data reduction, statistical analysis and to 
generate data displays. No attempt is made to analyze the scientific 
content of experiment results. The illustrations and plots herein 
are intended only to be graphical adjuncts to the textual description 
of the system. 

The authors have paraphrased from several unpublished documents 
and memorandums received from research personnel during the develop- 
ment, implementation and production phases of the effort. 


2.0 SCOPE OF THE DATA PROCESSING SYSTEM 


The overall effort for the processing of the OGO-F data involved 
three phases. The first phase consisted of the design, development 
and verification of the computer programming system to reduce and 
convert the digitized data to the phenomena the experiment was designed 
to observe. Basically, the primary outputs of this phase were the 
peak intensity values of definable lines from each processed scan of 
the six experiment sensors. The program design to accomplish this 
task was under continuous review during the entire contract to ensure 
that the most scientifically correct and usable results were obtained. 

The program was expanded and revised as analyses of the reduced data 
suggested new or better processing techniques or new avenues of 
scientific interest to be pursued. 

The second phase of the effort was the processing of the digitized 
data with the programming system developed in phase I. Seven months 
of scientifically interesting data were acquired. Each acquisition 
was stored as a file on the digitized tapes. Approximately 2,000 of 
these files were processed. The outputs from each file were 35mm film 
of aspect corrected counts versus wavelength for each scan of each 
sensor, a printout of the counts, wavelengths, collimator position, 
aspect and time for each frame of data, and punched data cards containing 
the peak intensities and background levels of predefined lines. The 
punched card data were ultimately stored on magnetic tapes in chronological 
order and were the inputs for the third phase of the effort. 

Phase three consisted of statistically analyzing and summarizing 
the peak data from phase II. Various listings and plots were generated 
and ultimately summarized results for each of the predefined lines were 
stored on two magnetic tapes. 

The three basic phases overlapped in time to a great extent. As the 
outputs from the production runs were analyzed, modifications and 
additions to more completely and accurately process the data were made 
to the programming system. Similarly, as the initial production punched 
card output was updated to magnetic tape, the programs to reduce, analyze 
and summarize these data were being designed, coded and verified. 
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3.0 EXPERIMENT TECHNICAL DESCRIPTION 


The experimental study which will utilize measurements from 
the spectrophotometer instrumented for the OGO-F spacecraft has, as 
objective, a determination of spectral and temporal variations of 
photon flux levels (radiated from the entire disk) in the extreme 
ultraviolet portion of the solar spectrum. Experimental emphasis, 
reflected in certain instrumental characteristics, is primarily 
directed toward what might be called an Intensity monitoring function 
rather than a detailed spectroscopic examination. Experimental 
objectives then are, repetitive, long-term observations, by means of 
a calibrated, scanning spectrometer having modest spectral resolution 
of intensity levels of those radiations in the extreme ultraviolet 
which are of aeronomical interest. 

The experimental package contains six collimator-grating spectro- 
meters in a parallel array which provide for an overall spectral 
scan ranging from 160 A to 1600 A. Throughout this wavelength range, 
the spectrometers have an effective resolving power of the order of 
100. By choice of grating constant, each of the six spectrometers 
covers a different spectral sub-range such that the total range 
scanned is as indicated above. The six Intervals are: 

Grating (lines/mm) Scan Limits (A) 


3600 

160-416 

2160 

267-695 

1800 

320-832 

1500 

384-998 

1200 

480-1248 

900 

640-1664 


These wavelength intervals, with the exception of the extremes, overlap 
to a considerable extent and as is apparent from the tabulation, 
gratings 2160, 1800 arid 900 1/mm include the Interval from 267-1664 A 
and gratings 3600, 1500, and 1200 1/mm cover the interval from 160- 
1248 A providing two nearly redundant groupings. 

Each of the six spectrometers consists of a planar grating, an 
analyzing slit system (collimator) and a detector plus associated 
electronics. The entrance aperture of each spectrometer is sufficiently 
large to allow uniform illumination of each grating for all orientations 
of the optical axis within ±2° from ideal point'lng along the mean 
solar vector. In the case of ideal spacecraft attitude, each point 
of the six gratings is illuminated at an 86® angle of incidence by 
solar rays which are 32 arcminutes divergent. This angular divergence 
generally determines the effective instrumental resolving power and 
leads to a value of about 100 for radiation from the entire solar disk. 
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Under illumination, spectra are produced in the respective planes 
of dispersion of the six planar gratings. Wavelength scanning is accom- 
plished with a six-channel collimator in spatial registration with the 
six gratings. The collimator drive is actuated' by a stepping motor 
which produces a 12 degree angular scan in 513 discrete steps. Motion 
of the collimator is equivalent to a scan of diffraction angles by the 
axial, high-transmission direction of the collimating slit system. 

At each step of the scan, dispersed radiation of the appropriate 
wavelength is transmitted by the collimator while radiation appearing 
at angles outside the acceptance angle of the slit system is rejected. 

The finite angular range of acceptance of the collimator' contributes to 
line broadening and, in part, determines the limiting resolution of 
the system. 

A channel electron multiplier (Bendix Channeltron) is located behind 
the exit slit of each of the six collimator apertures so as to receive 
the diffracted radiation which is appropriately transmitted at each 
position of the collimator. The channel multipliers are formed of 2 mm 
I.D. activated glass tubing bent in a 300 degree circular arc of 2.5 
inch diameter. With an overall voltage gradient of 3 kilovolts, the 
multipliers have a charge gain of approximately 10® which easily allows 
counting of output pulses produced by single photoelectrons generated 
at the input photocathode. The output terminal of the multiplier is 
operated at. high positive potential and is a.c. coupled directly into an 
integrated circuit differential amplifier which is located within the area 
circumscribed by the channel multiplier itself. These output signal pulses 
are fed, through a gate, to an 18-stage counter from which the stored 
information is subsequently shifted to the spacecraft data handling system. 

Data generated within the experiment which are of direct physical 
interest are (1) the angular location of the collimator relative to the 
diffraction gratings — this determines the wavelength for which the 
collimator transmission is maximum, and (2) the totalized count at a 
given collimator position for each detector — this determines the relative 
source intensity' (utilizing laboratory information) for various spectral 
locations. 

The collimator is step-driven (513 discrete steps per scan) through 
a 12 degree angular range (diffraction angles from 58® - 70®) and the 
step position is tabulated by a position counter which increments one 
unit for each driving pulse to the motor. The scan proceeds from step 
0 to step 513 at which point the direction is reversed to retrace from 
step 513 to step 0. Normally reversal of the drive direction is accom- 
plished by logic associated with the step-position counter, i.e. reversal 
is accomplished at steps 0 and 513. Mechanical switches are located at 
points which are about ten steps beyond the normal angular range to 
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provide reversing signals and counter reset information in the event 
of a severe lack of registration between the position counter and the 
actual mechanical scan position. In addition to the mechanical switches 
which provide end-of-scan reset information the electronic position 
counter is compared at points separated by 64 step positions with the 
output of an optical encoder (sunlight operated). The logic established 
by the encoder forces the counter into numerical coincidence with the 
unambiguous step position sensed by the encoder which is mechanically 
linked to the collimator drive system. Occasional corrections are 
necessary as the electronic counter may respond to certain spurious 
pulses (particularly those associated with the impulse commands) and 
fail to maintain registration with the actual mechanical scan position. 

The other measured quantities of interest are the photoelectron 
counting rates sensed by the detectors. At each of the 513 steps of the 
wavelength scan, the counter gates are cycled to sample the outputs 
from all six detectors, in sequence. In the detector electronics, each 
channel multiplier is followed by a pulse amplifier, the output of which 
terminates at a triple-input gate. Three amplifier outputs are switched, 
in sequence, to one eighteen-stage counter while the remaining three 
are switched to a second, identical coimter. The timing pattern for 
opening and closing the gates is determined by scan and sample logic 
within the experiment. Tor each of the 513 collimator positions, there 
will be six numbers which represent the accumulated count for the 
associated sampling interval, one number for each of the six spectrometers 
(one complete scan yields 3078 data points). Note that each step position 
gives intensity information at six different wavelengths. 

Solar Aspect System: 

The experimental package performs one additional measurement, that 
of experiment attitude (two-dimensional) relative to the mean solar 
vector. A bi-axial solar aspect system with an angular resolution of 
1/16 degree is mounted directly beneath the entrance aperture of the 
spectrometer. Aspect information is shifted to the spacecraft telemetry 
once every other main commutator frame. 

An aspect indicating unit is mounted on the sun-facing panel of the 
experiment package which will provide a coded readout of the orientation 
of the X- and z- axes of the package relative to the mean solar vector. 

The angular resolution of the device is 1/16 degree and the time resolu- 
tion, derived from alternate main commutator frames, is 288 msec at 8 
kilobits. Mechanically, the experiment is mounted with the optic axis 
parallel to the corresponding axis of the fine sun sensor aboard the 
spacecraft. This alignment is made with a tolerance of less than 1/4 
degree using optical techniques . 
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The aspect information is presented on the telemetry in sixteen 
bits. Of these sixteen bits, the first eight derive from the vertical 
reticle (parallel to the solar paddle x-axis) and the remaining eight from 
the horizontal reticle (parallel to the solar paddle z-axis) . Typical 
calibration data of binary readout (octal) as a function of angular 
displacement is given in the following table. 


Octal 


Angle (degrees) 


000 

001 

010 

020 

100 

177 

200 

300 

360 

370 

376 

377 


0.032 

0.096 

0.544 

1.056 

4.128 

8.160 

- 8.160 

- 4.064 

- 0.992 

- 0.480 

- 0.096 

- 0.032 
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4,0 FORMAT OF DIGITIZED DATA TAPES 


Experiment F-09 has been assigned four adjacent digital words in 
the spacecraft main frame, namely Words 49, 50, 51, 52, In addition 
the experiment utilizes two sub-commutator analog words to provide 
high voltage status information (Words 31 and 62), Two sequential 
main frames or eight digital words comprise the basic data unit of the 
experiment and alternate spacecraft frames are ordered in different 
data formats. Telemetry word and bit usage are as follows; 


- 

Word 49 


Housekeeping data 

Frame 1 < 

Word 50 


513-step position counter 

Format 1 

Words 51, 

52 

18-bit detector pulse 
counter — Counter No • 1 

Frame 2 j 

Words 49, 

50 

Solar Aspect Data 

Format 2 [ 

Words 51, 

52 

18-bit detector pulse 
counter — Counter No. 2. 


Frame (Format) 1, Word 49 

Bit No. 1 "1" Always positive — "one" 

2 "0" Counter identif ication-No. 1 


3 Scan Mode 

4 Drive Direction 

5 Mode Inversion 

6 Scan and 1 




7 Sample Modes J 

8 Detector 

9 Identification Code 


"0"-Long scan; "1" short scan 
"1"_0— >513; "0"-513~>0 
"1"-Real time; "0"-data storage 
j~ 00-Normal; 10-X4 step 
\ 01-X4 Sample; 11-X16 step 
01-Detector 1 

10- Detector 2 

11- Detector 3 

. 00-Stepping frame 
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Frame (Format) 1 - Word 50 

Collimator step position - A nine-bit binary code to identify steps 
0 to 512 (step 513 is also 0 following 512). 

Frame (Format) 1 - Words 51, 52 

Pulse Counter No. 1 

These two words provide telemetry for an eighteen-bit binary pulse 
counter which records pulses from Detectors No. 1, 2, and 3. 

Frame (Format) 2, Words 49, 50 

Bit No. 1 "1" 

2 " 1 " 

3 Aspect 

10 Sensor 

11 Aspect 
18 Sensor 

Frame (Format) 2, Words 51, 52 

Pulse Counter No. 2 

The two words provide telemetry for an eighteen-bit binary pulse 
counter which records pulses from Detectors No. 4, 5, and 6. 

In addition to the four main commutator words, three subcommutator 
words are stored in the basic frame structure. These are words 97, 98 
and 99. The subcommutators have a 128 step cycle. During each step 
each of the three words is monitoring a different spacecraft or experi- 
ment function. The current step (0-127) is encoded in word 65, one of 
the spacecraft ID words. Since each frame covers a time period of 288MS, 
a full 128 step subcommutator cycle is 18,43 seconds in the normal mode 
of operation. Consequently, each subcommutator parameter is updated 
every 18.43 seconds. The following seven subcommutator parameters 
were utilized in the analyses of the spectrometer data. 


Always positive - "one" 
Counter identification-No. 2 
Vertical Reticle - 
parallel to s/c X-axis 
Horizontal Reticle - 
parallel to s/c Z-axis 


4-2 


Subcom Word Code 


97-31 

HVl 

High Voltage Status 

for Detectors 1, 2, 3 

97-62 

HV2 

High Voltage Status 

for Detectors 4, 5, 6 

98-19 

D4 



98-20 

D5 



98-24 

All 



98-42 

AlO 



99-65 

E29 

General Thermistor Temperature 


In normal mode operation, the six detector outputs are sampled 
once per collimator step position and the complete sampling of outputs 
requires the utilization of eight successive telemetry frames. Since 
the experiment contains two binary pulse counters, two detector outputs 
may be sampled in parallel and the internal sampling program selects 
the pairs - Detectors 1 and 4; Detectors 2 and 5; and Detectors 3 and 6 ~ 
with the outputs 1, 2, and 3 switched to Counter No. 1 and outputs 4, 

5, and 6 to Counter No. 2. 

Consider eight sequential main commutator frames. Following Word 
52 in Frame 1, the gate from Detector 1 to Counter 1 is opened. This 
gate remains open until the beginning of Word 49, Frame 3 when the gate 
is closed. Data accumulated in Counter 1 during the above period are 
shifted out during Words 51, 52 of Frame 3 while Words 49, 50 contain 
housekeeping and collimator position data. At the end of Word 52, 

Frame 3, the gate is opened from Detector 2 to Counter No, 1 and is held 

open until Word 49, Frame 5 is reached, at which point accumulated data 
are shifted to the telemetry. The gate is then switched open for 
Detector 3 and is held open from the end of Word 52, Frame 5 until the 
beginning of Word 49, Frame 7 and these data are shifted during Words 49, 
50, 51, 52 of Frame 7. An analogous sequence is followed for Detectors 
4, 5, and 6 utilizing the intermediate frames, namely Frames 2, 4, 6, 
and 8. For example, the gate from Detector 4 to Counter 2 is open from 
the end of Word 52, Frame 2 until the beginning of Word 49, Frame 4, 
following which period, aspect system and Counter 2 data are shifted 
out to the telemetry. Following Word 52, Frame 7, all Counter No. 1 

data are inhibited and following Word 52, Frame 8, all Counter No. 2 

data are inhibited until Frame 9 (or 1), Word 52 and Frame 10 (or 2), 

Word 52, respectively, at which point the sequence begins over again. 

In the data inhibit period following Word 52, Frame 8, the collimator 
is advanced one step. For each detector, the gate is switched open for 
a duration of two frames less one 4-word shift-out interval, i.e. 288 - 
4.5 msec = 283.5 msec sample accumulation time. A collimator step 
occurs every eight frames or 1.152 sec and a long scan of 513 steps 
requires 9.85 minutes (a 65-step short scan requires 1.25 minutes). 

Figure 4.1 illustrates the basic data record. 
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DATA FORMAT FOR EXPERIMENTER 9 


CHARACTER 

REPRESENTATION 

1-2+42N 

D(33,j) 

Spacecraft Clock 

3-4+42N 

D(34,j) 


5-6+42N 

D(35,j) 


7-8+42N 

D(65,j) 


9-1 0+4 2N 

D(66,j) 

Spacecraft ID Words 

11-12+42N 

D(67,j) 


13-14+42N 

D(97,j) 


15-16+42N 

D(98,j) 

Subcoin Words 

17-18+42N 

D(99,j) 


19-20+42N 

D(107,j) 

SAI 

21-24+42N 

D(129,j) 

Fj^ Status field 


D(131J) 

F^ Status field 

25-26+42N 

D(49,j) 


27-28+42N 

D(50,j) 


29-30f42N 

D(5i,j) 

Main Comm Data 

31-32+42N 

D(52,j) 



33-34+42N 

D(132,j) 

Day of Year 

35-36+42N 

D(133,j) 


37-38+42N . 

D(134,j) 

Milliseconds of Day 

39-40+42N 

D(135,j) 


41-42+42N 

Spare 



0<N<63 Record #1 


64^N<127 Record #2 


Figure 4,1 
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5.0 DESCRIPTION OF COMPUTER PROGRAMS AND OPERATING INSTRUCTIONS 


Seven computer programs comprise the OGO-F Data Processing 
System, This section will contain a functional description of each 
and will describe the procedures for executing each program. 

5.1 Progranis RDLP and PKGROUP 

RDLP and PKGROUP are independent programs. However, they are 
executed as a two program system. RDLP, among other functions, 
outputs the peak intensities from each scan of data. PKGROUP analyzes 
the peaks and attempts to group the peaks from the multiple scans 
which represent the same line in the spectral range of each sensor. 

5.1.1 Functional Description 

Program RDLP contains several processing options which are defined 
via data cards. Among these options are multiple file processing, 
definition of parameters used in aspect and backlash correction, file 
selection, time period within a file to be processed selection, plot 
medium selection (35mm film or paper), scaling of the ordinate selection, 
listing of aspect corrected or uncprrected sensor counts and creation 
of a file for the interactive graphics program. The procedures for 
option selection will be detailed at the end of this section. 

The magnetic tape containing the digitized data is searched until 
the requested file is encountered. File identification information 
is extracted from the file header record. This information is printed 
as the first page of the listing and also is written as the first plot 
frame. Figure 5.1 is an example of the plot frame containing the file 
identification data. 

The data are unpacked one frame at a time until the millisecond 
of the day time field exceeds the requested start of processing time. 

The first frame of an eight frame message is located. This frame is 
identified by the detector index being equal to 0. The next successive 
seven frames are added to this first frame to comprise a message. A 
message contains all data acquired at one collimator position of a scan. 
The type of scan (long or short), scan direction, collimator position 
and time are noted. The six 18-bit sensor counts are extracted and 
are converted to count s/second. The three vertical and three hori- 
zontal aspect readings are extracted and converted to degrees using 
a calibration table. The vertical and horizontal aspect angles used 
in the aspect correction are weighted averages from the three vertical 
and horizontal components from each of three consecutive messages. 

The aspect angles for message N are calculated from the components in 
messages N-1, N and N+1. The components from message N are weighted 
more heavily than those from messages N-1 and N+1. If any components 
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contain fill data> the fill data are replaced by the average of 
legitimate components from the same message. Fill data can result 
from telemetry transmission problems or from problems during the 
digitization process. 

The scan position is corrected for vertical angle variations 
and mechanical backlash. Three corrected scan positions are calcu- 
lated for each message, and one is assigned to each pair of sensor 
readings with which it is time coincident. Each of the six sensor 
readings is corrected for horizontal and vertical aspect. The aspect 
correction is bypassed for any of the sensor readings which contain 
fill data. 

The wavelength for each sensor reading is calculated from the 
corrected grating position and the appropriate grating constant. If 
a graphics application file is being created, the time, grating position, 
six wavelengths, six uncorrected sensor values, six corrected sensor 
values, vertical and horizontal aspect and the three horizontal and 
vertical aspect components are written as one 29 word record. 

This message by message processing continues until a new long 
scan is encountered, six short scans have been accumulated or both 
power supplies are determined to be off. As each message is con- 
structed, the subcommutator position is interrogated. If the position 
is one of the seven of interest, the value in the appropriate sub- 
commutator word is converted to the desired units and saved. 

The scan start and stop time, and direction are saved. The 
corrected sensor readings are interrogated for eclipse. If any one 
sensor reading is negative, eclipse is assumed. If 50% or more of 
the sensor readings are between 5 and. 100 counts, eclipse is assumed. 

If the scan is in eclipse, the power supplies are off or scan contains 
fewer than 48 messages, the peak identification logic is bypassed. 

The aspect corrected sensor values for each sensor are now searched 
for peak values. A peak is defined by a minimum of three successive 
increasing intensities followed by a minimum of two intensities less 
than the highest intensity. Fill data are eliminated from this logic. 
The peak values, wavelengths and times are saved in arrays. - The peak 
data arrays are sorted so that all peaks are stored by increasing 
wavelength. 

The background intensity is calculated for each sensor in the scan. 
For short scans, the background is defined as the average or the minimum 
positive intensity and its adjacent intensities. For long scans, each 
sensor has unique criteria based on the wavelength spectrum covered by 
the sensor and the experiment configuration for that sensor. ' Following 
is a summary of these criteria. 
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IC 3600 


(1) Locate the minimum count between 370A-390A, 

(2) Average minimum with preceding and following counts, 
ignoring fill data. 

(3) Subtract minimum count from all peak values in the scan. 

IC 2160 

(1) Locate the minimxim count between 645A-665A.. 

(2) Average minimum with preceding and following counts, 
ignoring fill data. 

(3) Subtract minimum count from all peak values in the scan. 

IC 1800 

(1) Locate the minimum count between 645A-665A. 

(2) Average minimum with preceding and following counts, 
ignoring fill data. 

(3) Locate wavelength values and associated cotmts closest 
to 7 45 A and 81GA. 

(4) Solve exponential decay formula, I=Ae”^^, for A and X using 
two points from (3) . 

(5) From peaks with wavelengths ^ wavelength located in (1), 
subtract minimum count calculated in (2) . 

(6) For peaks with wavelengths ^ wavelength located in (1) , 
calculate I,' the background count, from I=Ae“^^ where x 
is the wavelength associated with the peak. Subtract I 
from peak count. 

IC 1500 

(1) Locate minimum count between 645A-665A. 

(2) Average minimum with preceding and following counts, 
ignoring fill data. 

(3) Locate wavelength values and associated counts closest 
to 810A and 880A. 

(4) Solve exponential decay formulas I"Ae“^^, for A and 1 
using two points from (3) . 

(5) From peaks with wavelengths ^ wavelength located in (1) 
and from peaks with wavelength ^ 915A, subtract minimum 
count calculated in (2). 

(6) For remaining peaks, calculate I, the backgroxmd count, 
from I=Ae“^^ where x is the wavelength associated with 
the peak. Subtract I from the peak count. 
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IC 1200 


(1) Locate mihimum count between 640A-660A, 

(2) Average minimum with preceding and following counts, 
ignoring fill data. 

(3) Locate wavelength values and associated coionts 
closest to 810A and 880A. 

(4) Solve exponential decay formula, I=Ae*^^, for A and ^ 
using two points from (3) . 

(5) From peaks with wavelengths _< wavelength located in 
(1) and from peaks with wavelengths >915A, subtract 
minimum count calculated in (2). 

(6) For remaining peaks, calculate I, the background . 
count, from I=Ae“^^, where x is the wavelength 
associated with the peak. Subtract I from the peak 
count. 

IC 0900 


. No background count is calculated for this sensor. 

The peak intensities in the scan for each sensor are corrected 
for background. The wavelength of the peak intensity determines 
whether the constant or exponential background, is subtracted for 
those sensors having both types of background. The peak data, 
backgrounds, times and wavelengths are written to a disk. The peak 
data are printed at the end of the listing and also are input to 
program PKGROUP. 

The time, scan position, corrected scan positions, corrected 
intensities, wavelengths and aspect angles are printed, one line 
per message. At the end of the printout for a short scan, the 
summation of all intensity values in the scan for each sensor is 
printed. For long scans, the summation is performed for each lOOA 
increment covered by the sensor. 

After the intensity summations for the last scan in a file 
have been printed, the subcommutator parameters accumulated during 
the processing of the file are listed. The parameters from one 128 
frame cycle are printed on one line. The time listed is from the 
frame containing the D4 parameter, which is for first parameter saved 
from each cycle. 

The subcommutator listing is followed by the peak intensity 
printout. For a file containing long scan data, the peaks selected 
from one scan are listed on two pages. The peaks are listed 
in ascending wavelength for each sensor. Peak intensities which 
had fill data in adjacent messages are noted with asterisks. 

The time of the peak to a thousandth of an hour, the peak wavelength. 
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intensity and intensity after background has been subtracted are 
printed. For short scan data, the printout differs in that peak 
intensities from six scans are printed. A scan number and background 
are printed for each short scan peak also. 

The plotted data consists of the six spectra from one scan per 
frame for long scan files and from six scans for short scan files. 

For long scan data, the 3600, 2160 and 1800 lines/mm sensors are 
all plotted with scaling from lOOA to 900A. The 1500 ^ 1200 and 
900 lines/mm sensors are all scaled from 400A to 1800A. This 
scaling facilitates the comparison or overlapping segments of the 
spectra. The ordinate scaling of the intensity values is a log 
base 10 representation with a range of 10^ to 10^ or 10^ to 10^. 

The abscissa scaling for short scan data is determined from the 
minimum and maximum wavelengths in the six scans being plotted, 
and consequently varies slightly from frame to frame. 

The constant and exponential background lines for each section 
of each spectrum are also plotted on the long scan frames. The 
horizontal and vertical aspect angles are plotted atop the 3600 
lines/mm sensor on a - 2 ° to +2® linear scale for both long and short 
scans. The ordinate range corresponds to the highest log decade 
from the intensity scaling. The time in hours, minutes and seconds 
is printed below the abscissa for the first and last data samples 
in the frame. Figure 5o2 is an example of a long scan plot frame. 
Figure 5.3 is an example of a short scan plot frame. 

At the conclusion of an execution of program RDLP, one or 
more files of peak intensity data have been created and stored on 
disk. One file exists for each file processed from the OGO-F 
digitized data tape. These files are processed by program PKGROUP 
which normally is executed at the termination of program RDLP. 

This dual execution eliminates the necessity of saving the file 
output from RDLP. The control card requirements for this dual 
execution are explained at the end of this section. 

The function of program PKGROUP is to group peak data 
representing the same predefined lines from the several spectral 
repetitive scans in a file. The grouping criteria are based on the 
closeness of the associated wavelength of a peak to predefined 
wavelengths and on the closeness of the peak intensity relative to 
the intensity of a major line from the same scan as the peak. Peaks 
which do not meet minimum standards are eliminated from further 
analyses. Peak intensities which vary significantly from the 
peaks representing the same line from other scans are statistically 
eliminated also. The relative intensity criterion does not apply 
to peaks from short scan files since in most cases the part of the 
spectrum investigated in a short scan does not include the major 
lines found in the complete sensor scan traversal. 
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The researcher has defined between 30 and 40 lines of interest 
for each of five of the six sensors. The 900 lines/mm sensor 
data does not undergo the peak grouping analyses. Figure 5*4 
lists the wavelengths and relative intensities for each sensor. 

The relative intensities are based on the ratio of the intensity 
of a given line to the intensity of one of the high amplitude lines 
in the part of the spectrum analyzed by the sensor. The following 
table lists the wavelength of the peak used for each sensor. 

Sensor Wavelength 


3600 

3 04 A* 

2160 

3 02 A* 

1800 

6 08 A 

1500 

6 08 A 

1200 

6 08 A 


* Both major peaks represent the same line. The difference 
is due to wavelength calculation inadequacies. 

These wavelengths and relative intensities are stored in arrays 
in program PKGROUP. The long scan data and short scan data are 
analyzed differently, so two separate descriptions will be given. 

A typical long scan file contains eight to ten scans. For 
each defined peak, the time, scan number, peak intensity, and peak 
intensity after background subtraction are read into arrays for 
each sensor. Consequently, all peaks in the file from one sensor 

are stored in one large array. Once all the peak data for the 

file has been stored, the analysis is performed on a sensor by 
sensor basis. Within a sensor, each scan undergoes the following 
processing. 

If the scan is a backward scan, 2A are added to each predefined 
wavelength and to the predefined major peak. A search is made for 

the major peak in the scan data. A peak within +5 A of the pre- 

defined major peak must be found, or the data from the sensor 
for the present scan is not included in the analysis. If the major 
peak is found, its intensity after background subtraction is saved 
to calculate the relative intensities of the other peaks in the 
scan. 

The wavelength of each peak in the scan is compared to the 
wavelengths of the predefined peaks for the sensor. If the 
wavelength of the peak is within +3A (3600 sensor) or H^A (2160, 
1800, 1500, 1200 sensors) of a predefined wavelength, the peak 
parameters are added to the group for that predefined wavelength. 
Because of the closeness of some of the predefined peaks, a peak 
from the data may be assigned to more than one predefined group 
using this gross criteria. Also, more than one peak from a 
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scan may be assigned to the same predefined line. To accommodate 
both situations, a measurement of closeness to a predefined peak is 
assigned each peak from the data. This measurement takes into account 
both the wavelength and relative intensity values of the peak from 
the data. The relative intensity is weighted roughly four times 
greater than the wavelength criteria. If more than one peak from a 
scan is assigned to the same predefined wavelength, the one having the 
largest measurement of closeness is tentatively assigned to the pre- 
defined wavelength. If the same peak is assigned to two predefined 
wavelengths and, if in both cases this peak has the greatest measurement 
of closeness to both predefined wavelengths from all peaks in its scan, 
the peak is assigned to that predefined wavelength for which its 
measurement of closeness. is larger. The peak is removed from the 
second predefined wavelength. If any other peaks from the same scan had 
also been assigned to the second predefined wavelength, the one with the 
largest measurement of closeness is assigned ,to that wavelength. This 
scheme eliminates the possibility of the same peak being assigned to two 
predefined wavelengths while ensuring that a given peak is assigned to 
the predefined wavelength to which it is best suited. based on the 
wavelength and relative peak intensity criteria. 

There also will be peaks which simply do not meet the ±3A or 
±4A criteria for any predefined wavelengths. These peaks are 
eliminated from any analysis and are listed as a group at the end of 
the printout. 

When all peak data from all scans in the file for a given sensor 
have been appropriately analyzed and uniquely assigned to the pre- 
defined wavelengths, the peaks in each predefined wavelength group 
undergo additional statistical analyses. The group intensity mean is 
calculated# Any peaks more than 1.1 standard deviations from the group 
intensity mean and also more than 35% from the mean are eliminated. 

Short scan peak data undergo similar analyses. However, rela- 
tive intensity is not calculated for short scan data. Consequently, 
the measurement of closeness criterion is based solely on the 
closeness of the data peak wavelength to the predefined wavelength. 

Also, the allowable variation from the group intensity mean is 
increased to 1.2 standard deviations. If a given predefined wave- 
length has fewer than eight peaks assigned to it from the entire 
file, the 35% from the group mean criteria is also used. 

For each predefined wavelength, the wavelength, relative 
intensity, group intensity mean before any peaks are rejected, 
group intensity mean after rejection and the standard deviation 
are printed. This summary is followed by the time, calculated 
wavelength, peak intensity, calculated relative intensity, scan 
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direction and reason for rejection if applicable for each peak 
assigned to the predefined wavelength. 

' For each file of data processed by PKGROUP, a BQ) card 
image file is created. The first card image contains file 
identification data. One or more cards images are written for 
each predefined wavelength. The exact number is determined by 
the number of peaks assigned to the wavelength. For each peak, 
the time, peak intensity after background subtraction, relative 
intensity and peak intensity before background subtraction are 
written. For short scan data, the relative intensity is always 0. 
These card image files are stored on magnetic tape in chronological 
order. 

5.1.2 Operating Instructions 

Programs RDLP and PKGROUP are independent programs. However, 
if peak grouping is desired, the two programs can be executed in 
the same computer run. RDLP is executed first and creates the data 
file read by PKGROUP. This dual execution eliminates the necessity 
of saving the intermediate data file created by RDLP. If RDLP 
processes multiple files, PKGROUP processes the data from all Of 
the files. This explanation assumes a basic understanding of the 
AFCRL CDC 6600 system control cards. Complete descriptions of 
these cards can be found in the AFCRL CDC 6600 users’ guide. 

The structure of the complete deck for compilation and 
execution consists of five sections. Each section is terminated by 
an EOR card (7/8/9 punched in card column 1). The first section 
contains the system control cards which perform functions such 
as compiling and executing the programs, requesting specific 
magnetic tapes and defining computer core and time requirements. 

The second section consists of the RDLP Fortran card deck. The 
third consists of data cards read by program RDLP. These cards 
contain program options and user definable variables used in aspect 
and backlash correction. The fourth section is the PKGROUP Fortran 
card deck. The fifth section is the data cards read by program 
PKGROUP. These are the predefined wavelengths and relative 
intensities. 


5. 1.2.1 System Control Cards 

Before using the AFCRL computer center CDC 6600, it is 
necessary to obtain an account number. Any users running with 
this number must be in the systems users' name library or attempted 
runs will abort because of the unrecognizable name on the job card. 
Also be sure that any magnetic tapes to be processed are submitted 
with the program deck. Refer to Figure 5.5 for this discussion. 

The cards in the system control cards section are bracketed by I. 


5-8 



Pigure 5.6 shows the two cards which must be submitted in the 
front of the punched card deck. The 6600 input card will be as 
shown. The only differences would be the input tape number and an 
increase of 250 seconds for each additional file processed beyond 
the first. The plot card is also standard. The estimated number 
of frames should be increased by 30 for each additional file 
processed in a multi-file run. 

Card #3, the blue job ID card, also requires an increase in 
the time field (T250) by 250 seconds for each additional file. 

The account number (7276) will be the one assigned by the computer 
center. 

Card //5 need be modified only by replacing the XXXX portion 
of (OGFXXXX) with the appropriate tape number. Note that the 
system control cards are terminated with an EOR card. 

5. 1.2. 2 RDLP Fortran Card Deck 

The RDLP Fortran Card Deck (II on Figure 5.5) begins with the 
program card for the main routine, RDLPY, and terminates with the 
END card from the last subroutine, LBCDI. This last END card must 
be followed by an EOR card. 

5. 1.2. 3 RDLP Data Cards 

Twenty data cards are required for processing one file of 
digitized data. Each additional file processed during one program 
execution requires two additional data cards. This section is 
III on Figure 5.5. 


Card 

Card 





// 

Columns 

Field 

Format 

Description 


1 

1-80 

Plot ID 

8A10 

A free formatted field. The data 

ap- 





pears at the beginning and end of 

the 





plots. User Name, laboratory and 
account number should be included t 


2 

1-5 

No. of Files to 
be processed 

15 

Low order digit must be in Column 

5. 

3 

2, 

Variable Number 

11 

Always 1 




K1 





3-31 

8 In Radians 
0 

F29.18 

Currently 1.00077 


4 

2 

Variable Number 

11 ■ 

Always 2 




K2 . 






6 Radians/Step 

F29.18 

Currently 0.00040771 


5 

2 

Variable Number 

11 

Always 3 




K3 





3-31 

a in Radians 
o 

F29.18 

Currently 1.501 


6 

2 

Variable Number 

11 

Always 4: 




K4 





3-31 

Vertical Zero 

F29.18 

Currently 0.0 


7 

2 

Variable Number 

11 

Always 5 




K5 





3-31 

Radian Conversion 

F29.18 

Currently 57.2958 


8 

2 

Variable Number 

11 

Always 6 




K6 





3-31 

Backlash Factor 

F29.18 

Currently 12.7 


9 

2 

Variable Number 

11 

Always 7 




K7 





3-31 

Scan Position 






Correction 

F29.18 

Currently 0.0 


10 

2 

Variable Number 
Vertical Aspect 

11 

Always 8 



3-31 

Correction K8 

F29.18 

Currently 4.0 


11 

2 

Variable Number 
Vertical Aspect 

11 

Always 9 



3-31 

Correction K9 

F29.18 

Currently 1.0 


12 

1-2 

Variable Number 
Vertical Aspect 

12 

Always 10 



3-31 

Correction KIO 

F29.18 

Currently 4.0 
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Card 

Card 





Columns 

Field 

Format 

Description 

13 

1-2 

Variable Number 
Vertical Aspect 

12 

Always 11 


3-31 

Correction Kll 

F29.18 

Currently 2.0 

14 

1-2 

Variable Number 
Vertical Aspect 

12 

Always 12 


3-31 

Correction K12 

F29.18 

Currently 1.0 

15 

1-2 

Variable Number 
Horizontal Aspect 

12 

Always 13 


3-31 

Correction K13 

F29.18 

Currently 0.25 

16 

1-2 

Variable Number 
Horizontal Aspect 

12 

Always 14 


3-31 

Correction K14 

F29.18 

Currently 0.0 

17 

1-2 

Variable Number 
Vertical Aspect 

12 

Always 15 


3-31 

Correction K15 

F29.18 

Currently 0.3279 

18 

1-2 

End of Variables 
Switch 

12 

Always 21 

19 

1-30 

Plot ID Data 

5A6 

File Description for Plots 
Ex; "Tape 0GF2OO File 4" 

20 

1-2 

File # to be . 
processed 

12 

Position of file on Input Tape. 


5-14 

Listing Start Time 

F10.3 

Time in seconds of day. Standard 
listing will not begin until this 
time is exceeded. Set to 0.0 usually 


17-26 

Listing Stop Time 

F10.3 

Time in seconds of day. Standard 
listing terminates when this time 
exceeded. Usually set to 99999.0 


29-38 

Plotting Start Time 

F10.3 

Same as 5^14 


41-50 

Plotting Stop Time 

F10.3 

Same as 17-26 


54 

Plot Medium 

11 

2 for Paper Plot 

3 for 35MM Film Plot 


55-59 Numerical Portion 

of Input Tape Label 15 


If Tape No. is 0GF196O, punch 1960 
in Cols. 56-59 
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Card 

// 


Card 


Columns 

Field 

Format 

Description 

62-64 

Log Scaling Option 

13 

10 for lOl - lO** 
100 for 102 - 10^ 

66-67 

Uncorrected Data 
Option 

F2.0 

1, for listing uncorrected 
sensor data. 0. otherwise. 

69-70 

Peak Determina- 
tion Option 

F2.0 

1, for peak determination. 
0, otherwise. 

72 

Graphics File 
Option 

11 

1 if graphics file is to be 
created. Blank otherwise. 
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Data cards 3 through 17 contain 15 variables used in various con- 
version and correction routines. Following is a summation of the 
mathematical formulas in which they are used. They will be referenced 
as K1 through K15. 

Variables Kl through K7 are used in the routine to correct the scan 
position (SCP) for angle variations and mechanical backlash. Three 
corrected scan positions (CWP) are generated for each SCP. The three 
CWP's are for sensors 1800 + 3600, 900 + 1500 and 2160 + 1200. The 
formula is: 

VERG 

CWP^ = ~ (ARCSIN(SIN(K3)-SIN(K3 k5^ ^4) + 

SIN(K1 + K2(.SCP + K6C1-CMK)))-K1) + K7 


1 , 2 , 3 

Vertical Aspect Angle for Sensor Pair X 

0 for Backward Scans 

1 for Forward Scans 

So in Radians 
5 Radians/Step 
“o Radians 
Vertical Zero 

Angular Degrees to Radians Conversion Constant 
Mechanical Backlash Correction in Scan Positions 
Arbitrary 

Variables K8 through K12 were used in the preliminary vertical aspect 
correction coefficient routine. However, they became superfluous when the 
routine was modified. All final production runs were made with the 
modified logic. The modified version uses only one of the variables read 
from data cards. The coefficient formula is: 

VERC^ = 1.0/ (1.0 + K15(VERG^)) 


where, 

X = 
VERG = 

X 

CMK = 


Kl = 
K2 = 
K3 = 
K4 = 
K5 = 
K6 = 
K7 = 


where. 


X = 1, 2, 3 

VERC^ = Vertical Aspect Correction Coefficient for Sensor Pair X, 
VERG^ = Vertical Aspect Angle for Sensor Pair X, 

K15 = 0.3279 

Variables K13 and K14 were used in the Horizontal Aspect Correction 
Coefficient routine. The horizontal correction coefficient formula is; 
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+ K14 


where, 

X = 1, 2, 3 

HORC^ = Horizontal Aspect Correction Coefficient for Sensor Pair X, 
HORG^ = Horizontal Aspect Angle for Sensor Pair X, 

K13 = 0.25 

K14 = 0.0. 

Cards 19 and 20 are repeated for each additional file processed 
during this execution. The last data card must be followed by an EOR 
card (7/8/9 multi-punched in column 1). 

5. 1.2. 4 PKGROUP Fortran Card Deck 

The PKGROUP Fortran Card Deck (IV on Figure 5.5) begins with the 
program card for the main routine, PKGRP, and terminates with the END 
card from the last subroutine, GROUP. This last END card must be followed 
by an EOR card. 

5. 1.2. 5 PKGROUP Data Cards 

The data cards for program PKGROUP consist of one card for each 
predefined wavelength as defined by the researcher. Although the number 
of predefined wavelengths varies for each sensor, a standard 40 cards 
are anticipated for each of the five sensors. Blank cards must be 
placed after the last legitimate wavelength card for any sensor with 
fewer than 40 predefined peaks to bring the total for that sensor to 
40 cards. The following table shows the number of blank cards needed 
for each sensor and the correct ordering of the cards. 


Sensor 

Cards 

Predefined Wavelengths 

Blank Cards 

3600 

1-40 

40 

0 

2160 

41-80 

30 

10 

1800 

81-120 

33 

7 

1500 

121-160 

33 

7 

1200 

161-200 

33 

7 


Within a sensor group, the cards must be in ascending order by wave- 
length. Columns 1-9 of each card contain the relative intensity in the 
form xxx.xxxxx. Columns 10-13 contain the corresponding wavelength 
number which must have its most significant digit starting in column 10. 
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The 200th card must be followed by an EOR card and an EOF card (6, 7, 8, 9 
multi-punched in column !)• 

5.1.2.6 Output From Programs RDLP And PKGROUP 

The printout from a dual execution contains for each file of data 
the file identification data, the variables from the data card input to 
RDLP, the scan by scan listing cf the sensor data, the subcoramutator data 
and the peak listing. Following the above printout for the last file 
processed is the listing of peaks grouped by predefined wavelength for 
each file processed. 

The plotted data consists of the file identification frame followed 
by one frame per long scan or by one frame per six short scans. Each 
frame contains a plot for all six sensors. For a multi-file processing 
execution, the plotted data for each file is preceded by a file identi- 
fication frame. 

The punched card data consists of a file header card followed by a 
variable number of cards containing the parameters for each peak grouped 
from the file. Multi-file processing executions will generate one of 
these decks for each file processed. However, the cards will be punched 
as one large deck with the file identification cards separating files. 

5.2 Background Application Program 

The punched card data from program PKGROUP were stored on magnetic 
tape. The approximately 3,000 files processed successfully were stored 
on nine tapes. Within a long scan file 169 variable length records 
were written. The size of the record was determined by the number of 
peaks from the file, grouped to a given predefined wavelength. Each 
short scan file contained 50 records since ten predefined lines accom- 
modated all possible ranges of the 16 types of short scans. The 
files were stored in chronological order. 

The following phases of analyses involved summarizing and plotting 
peak data for each predefined wavelength from multiple files. The 
parameter of interest to the researcher was the peak intensity after 
background subtraction. The criteria for background definition had been 
much more involved for long scan files than it had been for short scan 
files. For long scan files predefined wavelength ranges were analyzed 
for the best background definition. Some parts of the spectrum had 
background defined as an exponential function. For short scan files 
the background was defined simply as the minimum sensor intensity 
within each scan. Since both long scan and short scan data were to 
be part of the same analyses, it was necessary to standardize the 
background definition. The approach was to replace short scan backgrounds 
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with averaged long scan backgrounds from chronologically close files. 

For a given predefined wavelength in a short scan file, the 
background levels for the same sensor and wavelength from the five 
preceding and the five following long scan files were averaged. This 
average background replaced the preliminary background for all peaks 
at that wavelength in the short scan file. A new peak intensity after 
background subtraction was calculated for each peak in the short scan 
files. New magnetic tapes were created with the altered short scan 
data. The long scan data was not modified. All succeeding plot and 
statistical analyses programs used these newly created magnetic tapes. 

5.3 Peak Intensity List Program 


The Peak Intensity List Program prints the sensor, wavelength, 
average background and time, peak intensity before and after background 
subtraction and relative intensity for each peak for control card 
definable sections of one or more of the magnetic tapes containing 
peak data. Data selection options consist of a start and stop Julian 
day and hour, long and/or short scan files, individual or all five 
sensors and number of tapes to process. The deck setup for execution 
is as follows: 


Card # Column 


Description 


1 


2 . 

3 

4 

5 * 

6 

7 

8 


Orange or yellow 6600 input card. 

Core Memory = 50K 
C.P. = 300 seconds/tape 

P.P. = 300 seconds/tape 

Tape 1 = First Tape 

Tape 2 = Second Tape, if applicable 

DELOl, CM50000, T300, TP2. 7276 Delaney 
FTN. 


Request, Tape 1, HY. (Tape #/Noring/Delaney) 
Request, Tape 2, HY. (Tape ///Noring/Delaney) 
LGO. 

EOR Card (7/8/9 in Column 1) 


Fortran Card Deck 


EOR Card 


*If processing two tapes 
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Data 

Cards Column Description 


1 1-2 


3-4 

5-6 


10-16 

20-26 

2 1-9 

: 10-13 

201 

202 * 2-8 

12-18 

203 . 1 

204 1 


-1 = Process all Sensors 

1 = 3600 Sensor 

2 = 2160 Sensor 

3 = 1800 Sensor 

4 = 1500 Sensor 

5 = 1200 Sensor 

Number of tapes to be processed 
-1 = Process Long and Short Scans 

0 = Process Short Scans only 

1 = Process Long Scans only 

Start time for processing in Julian 
day, hour and hundredths of hour 
Ex: JUN 29 at 3«70 hours = 1790370 

Stop time for processing in Julian 
day, hour and hundredths of hour 
Relative Intensity T For predefined 
Wavelength \ Wavelengths 


Start time for processing for second 
tape if applicable 

Stop time for processing for second 

tape if applicable 

7/8/9 

6/7/8Z9 


*0nly needed if second tape being processed 
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5.4 


Relative Intensity Analysis Program (OGOREL) 


A second parameter for each of the predefined wavelengths 
used for the peak intensity grouping analysis was the standard 
relative intensity of each predefined peak to one of the major 
peaks in the section of the spectrum covered by the sensor. 

The major peak for each sensor was chosen for its consistent 
identification based on amplitude and positioning within a scan 
sweep. A secondary major peak was also chosen in the event that 
the primary major peak was not identifiable because of noise or 
fill data. The relative, intensity for each predefined wave- 
length for a sensor was simply the ratio of the peak intensity 
of that line to the peak intensity of the major line. The 
standard relative intensity was determined by averaging the actual 
ratios from several files of data. This relative intensity 
criteria added a second dimension to the peak grouping logic which 
was based on the data from the scan in which a peak was recorded. 
Program OGOREL was designed to statistically compare calculated 
relative intensities from many long scale files to the standard 
relative intensity for each predefined wavelength. 

Program OGOREL utilizes the peak amplitude data from the 
peak grouping analyses. For long scan files, the calculated relative 
intensity is one of the parameters stored for each peak. OGOREL 
selects the data from magnetic tape from long scan files only 
since the relative intensity is not calculated from short scan 
data. A maximum of 100 long scan files can be accomodated. The 
relative intensities and times for all peaks for each predefined 
wavelength for all long scans are stored by sensor and predefined 
wavelength. Thus, if 50 long scan files are read, each with three 
peaks at a given predefined wavelength, a total of 150 relative 
intensities are used for that predefined wavelength in this 
analysis. 

The accumulated relative intensities for each line undergo 
the following processing. Any relative intensities which are 
less than one tenth or greater than fifty times the standard 
relative intensity for the line are eliminated. This gross 
filtering rejects data from scans where the major peaks are 
ill defined, and consequently the relative intensities from that 
scan are distorted. The estimated standard deviation is calculated 
for the remaining relative intensities by: 


where , 

ESD 

N 

X 

o 


ESD = 


T (X.-X 

1=1 X 0 

N-1 


Estimated Standard Deviation for the line 
Number of relative intensities 
Mean Relative Intensity 
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Any relative intensities outside +1*1 CESD) of the mean are 
rejected, and the mean is recalculated. The sensor, wavelength, 
standard relative intensity, number of samples, average relative 
intensity before any samples are rejected and average relative 
intensity after rejection are printed. Each sample and its time 
is then listed. Those eliminated by the gross filter are flagged 
with a dollar sign; those eliminated by the standard deviation 
criterion are flagged with an asterisk. At the end of the analyses, 
a summary is listed for each line for each sensor. The wavelength, 
standard relative intensity, calculated mean relative intensity 
after rejection, the % difference between the two and the number 
of samples in the analyses are printed for each predefined line. 

The deck setup for an execution of OGOREL is as follows: 


CARD DESCRIPTION 

1 Orange or yellow 6600 input card. 

Core Memory = 114K 

C.P. = 180 seconds 

P.P. = 500 seconds 

Pages = 100 

Tape 1 = Tape # for peak data 

2 DELOl, T180, CM114000, TPOl. 7276 Delaney 

3 FTN. 

4 REQUEST, TAPEl, HY. (TAPE # /NORING) 

5 SETCORE. 

6 LOAD, LGO. 

7 EXECUTE. 

8 7/8/9 IN COLUMN 1 

9 OGOREL FORTRAN CARD DECK 


. . 

7/8/9 IN COLUMN 1 

The 200 predefined peak and standard relative intensity 
intensity cards as described in Section 5.1 
7/8/9 IN COLUMN 1 
6/7/8Z9 IN COLUMN 1 

5.5 Peak Plot Program 


The peak plot program processes the peak data after the back- 
ground has been subtracted. Its purpose is to illustrate the peak 
amplitude profile for each predefined line over a variable time 
period. Both long scan and short scan files are utilized. A plot 
is generated for each of the 169 wavelengths. Five sets of plots 
are stacked vertically on one plot frame to reduce the plotter 
requirements. The time period selected is an input variable, but 
practically ranges from two to four day periods, ^ program 

execution, all abscissa scaling is the same. The ordinate scaling 
for each of the 169 lines is control card definable. Either 35mm 
film or paper plots can be requested via data card option. 
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The magnetic tape containing the requested peak data is 
searched until the first file within the desired time frame is 
encountered. The peak intensity after background subtraction and 
the time is extracted from all peaks and stored by predefined 
wavelength. The sensors are plotted in order by descending grating 
lines/MM. Within each sensor, the data are plotted by ascending 
wavelength, five per frame. The first wavelength of a sensor 
generates a new plot frame. 

Each frame contains sensor, wavelength, and time period identi- 
fication. A vertical line is drawn at the 0 hour of each day. Each 
of the five abscissas are labelled depending on the scaling chosen 
for each wavelength. Each peak value for a wavelength is plotted 
as an individual point. Short scan data can usually be differentiated 
from long scan data because of the high density of points clustered 
together. If a value exceeds the ordinate scaling limit, it is 
entered in the plot above the one for its wavelength with ordinate 
position determined by the scaling for its wavelength. It will not 
be confused with data from the upper plot since each successive plot 
alternates the two characters, "X" and "0". Consequently, if a 
particular plot is utilizing the symbol "X", all of its values in 
the plot above it will be "X"'s in a plot with predominately "0" 
symbols. Values in the fifth or uppermost plot of a frame which 
exceed that plot’s maximum are not plotted to avoid abnormal 
termination of the program execution. Figure 5.7 is an example of 
a peak plot frame. 

Following is the deck setup for the peak plot program. 

Card Column Description 

Orange or Yellow 6600 Input Card 
Core Memory = 77K 

C.P. = 300 seconds 

P.P. = 500 seconds 

Plot = 

Tape 2 = Tape // of peak data 

Tape 39 = CRTPLTS or PENPLOTS 

GREY PLOT CARD 

DELOl, T300, CM77000, TP2. 7276 Delaney 
FTN. 

SETCORE. 

ATTACH, CRT, NEWCRT, CY=3, MRal. OR, 
ATTACH, PEN, PENPLOTS, Mk:±l. 

REQUEST, TAPE2, HY. (PEAK TAPE ///NORING) 
REQUEST, TAPE39, HY. (CRTPLTS /DELANEY) OR, 
REQUEST, TAPE39, HY. (PENPLOTS/DELANEY) 
LOAD, CRT. OR, 

LOAD, PEN. 

LGO. 


2 

3 

4 

5 

6 

7 

8 

9 

10 
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Card 


Columns 


Description 


11 

12 


Data 

Cards 

1 


2 


3 


4 

• 

203 

204 

205 


1 7/8/9 

PEAJC PLOT FORTRAN CARD DECK 


7/8/9 

data CARDS 


f/8/9 

6/7/8/9 


Columns Description 

1-4 Total number of hours of data to be 

processed. 

EX; 4 Day Plot = 96 
2 Day Plot = 48 
6 Plot Medium Option 

1 = 35mm film 
. 2 = paper 

1-30 Free Formated Plot ID 

1 . 10 . 20 

EX: Delaney 7276 OGO-F 

1-7 Start Julian Day and Hour (JJJHHMM) 

N.B. The Julian day punched is 
actually one less than the desired 
day. 

EX; For January 30, 00 hours 
punch 029 0000 

10-16 Stop Julian day and hour 

N.B. The Julian day punched is 
actually one less than the desired 
day. 

The 200 predefined wavelength 
cards as described in section 

5.1. 

1 7/8/9 

1 6/7/8Z9 
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5.6 


Statistical Summary Program (SATAV) 


The function of program SATAV is to summarize and to statistically 
analyze the peak amplitude after background subtraction and relative 
intensity data. For each predefined wavelength, this analysis is 
performed for each requested Julian day. Within each day, analyses 
are performed on each three hour segment, each six hour segment, each 
twelve hour segment and the entire twenty-four hour period. For 
each time segment, the number of samples, peak amplitude mean and 
the unbiased estimate of standard deviation are calculated. For 
each time segment, the same quantities are calculated for the non- 
zero (long scan data only)’ relative intensity ratios. Any non-zero 
relative intensities within the time segment which are not within 
(1.1) X (unbiased standard deviation) are eliminated, and the mean 
relative intensity and unbiased estimate of standard deviation are 
recalculated from the remaining relative intensities. 

The output consists of a listing and a magnetic tape file. For 
each day processed, the following is printed for each of the 169 
predefined wavelengths. The sensor, wavelength in angstroms and the 
date of data acquisition are printed followed by the time (hour and 
hundredths of hour as an integer number), peak intensity, background 
and relative intensity ratio for all peaks in the entire twenty-four 
hour period for the predefined wavelength. For each of the eight 
three-hour segments, for each of the four six-hour segments, for 
each of the two twelve-hour segments and for the entire twenty-four 
hour period, the following are printed; the beginning and end hour 
of the time segment, the number of peak amplitude samples in the 
time segment, the mean peak amplitude (i^y)> the unbiased standard 
deviation (a_), + o_, - a_ the number of non-zero relative 

intensity ratios after elimination, the mean relative intensity ratio 
after elimination and the unbiased estimate of standard deviation 
after elimination. These same quantities for each time segment are 
written on the magnetic tape file, one record per Julian day for 
each of the 169 predefined wavelengths. Figure 5.8 is an example 
of the SATAV printout. 

Following is the deck setup for program SATAV. 



1 Orange or Yellow 6600 input card. 

Core Memory = 120K 

C.P. =30 seconds for each day 

to be processed. 

P.P. = 100 seconds for each day 

to be processed. 

Tape 1 = Tape # of peak data 

Tape 2 = Tape # of summary output 

2 DELOl, T30*, CM120000, TP2. 7276 Delaney 
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Card 


Description 


3 

4 

5 

6 

7 

8 
9 


FTN. 

REQUEST, TAPEl, HY. (PEAK TAPE ///NORING) 
REQUEST, TAPE2, HY. SUMMARY TAPE ///RINGIN) 
SETCORE. 

LGO. 

EOR CARD (7/8/9 in Column 1) 

SATAV FORTRAN CARD DECK 


EOR CARD 

DATA CARDS (explained below) 
EOR CARD 

EOF CARD (6/ 7/ 8/9 in Column 1) 
* Add 30 for each additional day processed. 


Data 

Cards Columns Description 


1 1-7 


9-15 

2 2-5 
7-10 

12-15 

17-20 

22-25 

3 10-13 

i ft 

202 10-13 


First day to be processed in Julian day, 
hour and hundredths of hour as an integer 
number. Format is: JJJHH.HH (Decimal point 

assumed) 

Last day to be processed in Julian day, hour 

and hxmdredths of hour as an integer number. 

Always = 3600 

Always = 2160 

Always = 1800 

Always = 1500 

Always = 1200 

Predefined wavelengths in angstroms. 

Same format as described in Section 5.1. 
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5.7 Statistical Siunmary Plot Program (SUMPLOT) 

The function of program SUMPLOT is to plot the 12 or 24 hour 
peak intensity averages versus time for each predefined wavelength 
over a data card defined time period. The plot frame format is 
similar to that described for the peak Plot Program in section 5.5 
Five sets of axes, one for each of five predefined wavelengths, 
are stacked vertically on one plot frame. Format differences 
include abscissa scaling which is ten days/inch for SUMPLOT and 
the addition of error bar lines for each sample plotted. , The 
error bars are the average peak intensity plus or minus the unbiased 
standard deviation for the samples in the 12 or 24 hour time segment. 
The ordinate scaling for each of the 169 predefined wavelengths is 
defined via data cards. 


This program was designed primarily to aid in the analysis of 
the statistically summarized data from the entire six month data 
base. Consequently, the abscissa scaling is a constant 10 days/ 
inch. The length of a plot frame is solely determined by the total 
number of days being plotted. For computer turnaround purposes, 
the program only plots summary data for one sensor during one execu- 
tion. Consequently, five runs are necessary, to obtain plots for 
all 169 predefined wavelengths for a desired time period. Figure 5.9 
is an example of a frame of plotted summary data. The scaling of 
this example was modified for publication purposes. 

The input data to program SUMPLOT is generated by program SATAV. 
Program SATAV created nine magnetic tapes for the entire data base. 
Since an execution of program SUMPLOT usually utilizes the entire 
span of acquisition, an intermediate step was taken to merge the 
nine tapes from SATAV onto two magnetic tapes. Consequently, 
program SUMPLOT has been designed to read both tapes. 


The deck setup for program SUMPLOT follows. 


Card 


Description 


1 


Orange 

Core 

C.P. 

P.P. 

Tape 

Tape 

Tape 

Plot 


or Yellow 6600 Input Card, 
memory = 77K 

= 500 seconds 
= 1000 seconds 

1 = Tape # of chronologically 

first summary tape. 

2 = Tape # of chronologically 

second summary tape. 

39 = PENPLT or CRTPLT 

= Check, box on request card 
if plot is required. 
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Card 


Description 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


GREY PLOT CARD. 

DELOl, CM77000, T500, TP2. 7276 DELANEY 

FTN. 

REQUEST, TAPEl, HY. (1st TAPE/NORING) 
REQUEST, TAPE2, HY. (2nd TAPE/NORING) 

SET CORE. 

ATTACH, PEN, NEW ON LINE PEN. OR, 

ATTACH, CRT, NEW ON LINE CRT. 

LOAD, PEN. OR, 

LOAD, CRT. 

LGO. 

EOR CARD (7/8/9 in Column 1) 

SUMPLOT FORTRAN CARD DECK 


EOR CARD 

DATA CARDS (DESCRIBED BELOW) 


EOR CARD 

EOF CARD (6/ 7/8/9 in Column 1) 


Data 

Cards Columns Description 


1 1-30 

2 1-3 
8-10 

16 

18-21 

3 9-12 
16-20 


30 Character Plot Identifier 
EX: Delaney 7276 OGO-F 

Start Julian day of year 
Stop Julian day of year 

1 = Paper Plot 

2 = 35mm Film Plot 

Sensor to be processed. Must be one of the 
following; 1200, 1500, 1800, 2160 or 3600. 
First predefined wavelength for selected 
sensor 

Maximum value for ordinate scaling for first 
wavelength 


42 9-12 

16-20 


*40th predefined wavelength for selected 
sensor 

Maximum value for ordinate scaling for 40th 
wavelength 
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Data 

Cards 


Columns 


Description 


A3 1 EOR 

A 4 1 EOF 


*If fewer than 40 predefined wavelengths for selected sensor, 
blank* cards must be added at end to bring total to 40, 
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FILE IDENTIFICATION DATA 


TRPE NO. QOOF 140 
FILE NO. 10 
EXPERINENT NO. s 9 
SATELLITE ID s 69511 
DATE - JUN 9.1969 (6/9/69) 

START TIAEs 6HRS. 36N1NS. S2SECS. (29912 SECONDS OF OAT 160) 

STOP TINE: 7HRS. 95A1NS. 81SECS. (27961 SECONDS OF OAT 160) 

0R6IT NO. s S3 
STATION NO. s 21 

PLOTTING FROn 0>0»0.00 (0.00 SECS.) TO 24»0»0.00 (06400.00 SECS.) 



Figure 5.1 
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LONG SCAN PLOT FRAME 









//y/ 



SHORT SCAN PLOT FRAME 



WAVELENGTH & RELATIVE INTENSITY VALUES 
3600 Lines/mm 2160 Lines/mm 1800 Lines/nm 


Predefined 

Wavelength 

Relative 

Intensity 

Predefined 

Wavelength 

Relative 

Intensity 

Predefined 

Wavelength 

Relative 

Intensity 

175 

.02456 

302 

1.00000 

348 

.06231 

177 

.02480 

314 

.03592 

353 

.06362 

180 

.03810 

320 

.03208 

361 

.16710 

185 

.01540 

335 

.19720 

368 

.12940 

188 

.03305 

348 

.06247 

377 

.01218 

193 

.03603 

361 

.12790 

391 

.02288 

195 

.04110 

368 

.10340 

403 

.02591 

202 

.03826 

377 

.01311 

417 

.04025 

211 

.03566 

391 

.01928 

423 

.01378 

219 

.02014 

403 

.02170 

430 

.01220 

227 

.02055 

417 

.03396 

436 

.02915 

234 

.01304 

423 

.01224 

464 

.06363 

243 

.02289 

430 

.01486 

500 

.21030 

252 

.02283 

436 

.02447 

512 

.07116 

256 

, .08187 

447 : 

.01575 

521 

.09601 

265 

.03127 

464 

.04621 

529 

.02553 

270 

.01716 

500 

.12120 

537 

.05078 

274 

.03521 

512 

.03132 

554 

.15930 

277 

.01321 

521 

.05426 

568 

.16880 

284 

.18950 

529 

.01122 

586 

.53220 

291 

.01363 

537 

.02830 

608 

1.00000 

295 

.01449 

554 

.07581 

630 

.46950 

304 

1.00000 

568 

.06664 

640 

.02524 

314 

.02004 

584 

.24110 

672 

.19260 

320 

.02102 

600 . 

.02591' 

686 

.05761 

335 

.17320 

608 

.34870 

695 

.06388 

342 

.01837 

630 

.19630 

705 

.12780 

348 

.04670 

672 

.05073 

722 

.12920 

353 

.03670 

688 

.07939 

737 

.10420 

361 

.11040 

710 

.04443 

770 

.18400 

368 

.10020 



790 

.18970 

377 

.00930 



825 

.03135 

381 

.00185 



838 

.21900 


391 .01591 
397 .00175 
403 .01502 
409 .00544 
414 .00369 
419 .02916 
427 .01086 


Figure 5.4 
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WAVELENGTH & RELATIVE INTENSITY VALUES 


(cont 'd) 


1500 Lines/nnn 

1200 Lines/mm 

Predefined 

Relative 

Predefined 

Relat ive 

Wavelength 

Intensity 

Wavelength 

Intensity 

417 

.05270 

537 

.05983 

430 

.01624 

554 

.18360 

436 

.03860 

568 

.16270 

445 

.01941 

584 

.63040 

. 464 

.07514 

608 

1.00000 

500 

.24120 

630 

.56240 

512 

.07874 

640 

.02468 

521 

.11540 

672 

.17630 

529 

.02274 

686 

.07154 

537 

.06133 

705 

.15500 

554 

.18400 

722 

.13360 

568 

.16780 

737 

.10040 

584 

.59990 

770 

.23450 

608 

1.00000 

790 

.25080 

630 

.55020 

834 

.28700 

640 

.02545 

853 

.16360 

672 

.17040 

875 

. 03814 

686 

.06943 

912 

.73140 

705 

.15070 

932 

.26800 

722 

.12480 

951 

.20310 

737 

.09886 

977 

1.22680 

770 

.23130 

998 

.37840 

790 

.25800 

1005 

.40700 

834 

.29880 

1011 

.41540 

853 

.15550 

1026 

.95780 

875 

.04533 

1057 

.17680 

912 

.73320 

1081 

.34040 

932 

.24000 

1097 

.18940 

951 

.20620 

1108 

.36550 

977 

1.21680 

1138 

.33980 

992 

.23110 

1170 

.81180 

1002 

.33160 

1254 

.38400 

1011 

.35340 

1266 

.66310 
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Card # 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


RDLP AND PKGROUP DECK SETUP 


Contents 


,—ORANGE OR YELLOW 6600 INPUT 
I GRAY PLOT CARD 
I DELOl, CM165000, T250, TP2. 
i FTN. 

i REQUEST, TAPE 3, Hi, S. 

I REQUEST, TAPE 39, HY. 

^ SETCORE. 

4 ATTACH, CRT, CRTPLOTS, MR^l 
“■ LOAD, CRT. 

LOAD, LGO. 

EXECUTE. 

FTN. 

RFL, 130000. 

LOAD, LGO. 

. EXECUTE . 

^7/8/9 


pRDLP FORTRAN DECK 

I : 

TL 


L 


7/8/9 

RDLP DATA CARDS 


L: 


7/8/9 

— PKGROUP FORTRAN. DECK 

'^ 11^19 

PKGROUP DATA CARDS 
3T' 7/8/9 
t-6/7/8/9 


CARD 

7276 DELANEY 

(OGFXXXX/NORING) 

(CRTPLT/RING/DELANEY) 


Figure 5.5 
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CDC 6600 INPUT CARD AND PLOT REQUEST CARD 


AFCRL COMPUTER CENTER SYS 

ITEM II 

€ 

»600 INPUT CARD ^ 

PROGRAMMER 

b>u.AKii-.y 

JOB ID 

O J 

TEL 

CLASS. 

PRODUCTION er 
CHECK OUT a 


STATE INPUT 
DISK PACK 
PT» PAPER TAPE: 
BLANK B MAG TAPE 

TAPE/DISK 
NAME (LA8EU 

LIST 

EST. NO. PAGES ^ 0 _ 

CORE MEMORY REQUIRED 1 

OCTAL 




CREATE FILES 
P.F. □ 

COMMON 1 1 

ESTIMATE TIME 
DECIMAL 

c p /i* / ) srr 

?, 7 






P. P. / 0 0 SEC, 




SPECIAL INSTRUCTIONS: 

PROGRAM 

DECK 

RETURNED 

(!□ 

TAPES 

COPY 1 J 

DUMP 1 . J 

VERIFY 1 ,] 

OUTPUT 

PUNCH 1. 

PLOT (, ^ 


AFCRL rOMM 60 lJUN 71) OtTV.I «-7a 
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STATISTICAL SUMMARY PROGRAM PRINTOUT (SATAV) 




6.0 OGO-F PRODUCTION MASTER FILE LOGGING SYSTEM 


The OGO-F satellite transmitted data for almost seven months. An 
average of 30-35 acquisitions were digitized for each universal day. 
Generally, two magnetic tapes contained acquisitions for one universal 
day, one acquisition per physical file. These data were stored on over 
500 magnetic tapes. Data for some universal days required a third magnetic 
tape for various reasons, and some magnetic tapes contained a few acquisi- 
tions from each of several universal days. Consequently, it was evident 
that information pertaining to the contents of the digitized data tapes 
and the types of processing accomplished with the data must be made avail- 
able to the research personnel and to the production personnel. The 
research personnel had to know if data were available in periods of 
scientific interest and if the data had been processed with the RDLP, 
PKGROUP, or the OGO Graphics programs. The production personnel had to 
know the contents of each digitized tape in order to proceed in a 
logical manner and to be able to make intelligent estimates of the com- 
pletion of processing. 

Generally, within a physical tape, data was stored chronologically, 
universal time being the only parameter unique to each acquisition. 
Therefore, a chronological record of the data was created. 

The following subsections explain the mechanics of the logging 
system from both a functional and practical viewpoint. 

6.1 Sources of Data 

6.1.1 OGO-F Data Tape Packing Sheet 

Accompanying each magnetic data tape received from NASA was a packing 
sheet. The packing sheets contained such information as the NASA assigned 
tape number, start and stop times, recording station, recording dates, 
and total number of files. Pertinent information was extracted from the 
individual records on the packing sheets and used as a basis for the 
production logging system. Figure 6,1 is an example of a packing sheet. 

6.1.2 OGO-F Data Parameters from Research Personnel 

Research personnel supplied us with data they wanted incorporated 
into the Master File. The parameters included experiment power supply 
on/off times, scan code, scan time, flare number, flare type, flare 
start and stop times, and flare maximum time. Additional data supplied 
by research personnel consisted of start and stop times for eclipse periods. 
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daily 10.7 CM flux averages and daily averages. 

6.2 Master File Data Record 

The Master file consisted of 100-word data records (see figure 6.2), 
one record for each data acquisition. Information extracted from the packing 
sheet included date of recording, start time, stop time, and file number. 

The tape number was identified by labeling the tape ”0G0Fxxxx”, where 
xxxx was a four digit number. The tapes were labeled in order by date of 
acquisitions in increments of 10. Tapes received after the effort which 
chronologically fit between the original tapes, were labeled accordingly, 
having numbers ending in a 3, 5, or 7. 

The parameters which the research personnel supplied were incorporated 
into the Master File data record according to the start time of the event. 

If the start time was found to be between two recording periods, the infor- 
mation was stored in the data record which preceded the start time of the 
event. For each 100 word data record, a punched card was generated for 
updating production information to the Master File. Figure 6.3 is the 
format for the punched card. 

When data processing was completed and verified, for a file, a file 
status code, date processed, film number, and a film code were added to 
the punched card for that file and updated to the data record for that 
file in the Master File. (See Section 6.4o) 

6.3 OGO-F Master File Printout 

The Master File printout is a computer listing generated to aid 
research personnel in determining data acquisitions of interest. Perti- 
nent information was extracted from the 100-word data record for each 
file which comprise the OGO-F Master File. The parameters include start 
and stop times of data acquisitions, tape and file numbers, daily 10.7CM 
flux averages, daily A averages, power supply on/off times, scan code, 
time of scan change, f?are event times, flare event types, and eclipse 
periods. Figures 6.4 and 6.4.1 are examples of the OGO-F Master File 
printout. Figure 6.5 is an explanation of each field on the Master 
File printout for 1969. 
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6,4 Updating Production Master File 

The production master file was updated periodically for the purpose 
of keeping accurate account of files processed with the RDLP and PKGROUP 
programs. 

Verification of production files consisted of reviewing the RDLP 
output listing, reviewing the output listing and punched card peak data 
from PKGROUP, and reviewing CRT film for both quality and correctness. 
When this task was completed, a file status, date processed, film number, 
if applicable, and film code, if applicable, were added to the punched 
card (see section 6.2) corresponding to that acquisition on the Master 
File. 


A number of data files could not be processed for various reasons. 

A numerical code (see Figure 6.7) was added to the file status field 
of the punched card which corresponded to these acquisitions on the 
Master File. Figure 6.6 is an example of the punched card format used 
for updating this function to the OGO-F Master File. 

The punched cards accumulated for one update period were sorted 
chronologically. For these acquisitions, the file status, date processed, 
film number, and film code (CRT/PEN) were added into the appropriate 
100-word record. Figures 6.8 and 6.8.1 are examples of the OGO-F 
production log status. Figure 6.9 is an explanation of each field on 
the production log printout. 

6.5 OGO-F Chronological Printout - 15 -Minute Intervals 

The computer listing consisted of data parameters extracted from the 
100-word record. Each day was printed in two 12-hour segments. Each 
12-hour segment is divided into 48 15-minute periods. Each 15-minute period 
was represented by one line of printout. Basically, any of the data, 
parameters, or times which occur, begin, or continue during any part of 
the 15-minute period are printed in a designated field for that particular 
line of printout for the 15-minute period. Change of power supply status 
or scan type is noted. Figure 6.10 is an example for the chronological 
printout in 15-minute intervals. 
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OGO-F DATA TAPE PACKING SHEET 



OGO EXPERIMENTER 

STA TAPE file RECORDED EDIT EDIT R FF START STOP DfCoH ExP 

YR MO OY TAPE file HR MN HR MN RUN 9 

T F 
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1 69 
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16 

1 78 1 

1 3 
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27 

691991 
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SKA 
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1 69 
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16 

1781 

I 3 
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□ 

27 

0 
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691991 
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SKA 

3H9 

1 69 
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16 

1781 

1 3 
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S*7 
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T8 
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69199 1 
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*) 
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1 69 
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|78r 
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0 

2 

10 
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19 

691991 

1 

5 

act 
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1 69 

7 

18 

1786 

1 3 

0 

2 

15 

3 

B 

69 1 99 1 


6 

ACT 

256 

1 69 

7 

16 

1786 

1 3 

0 

3 

50 

3 

57 

691991 


7 

ROS 

233 

1 69 

7 

18 

1788 

1 3 

0 

3 

57 

5 

31 

691991 

1 

B 

ROS 

233 

1 69 

7 

1 8 

1789 

1 3 

0 

5 

33 

5 

50 

691991 


9 

ROS 

2 33 

1 69 

7 

18 

1788 

1 3 
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59 
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22 

691991 
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ROS 

2 33 
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16 

1788 

1 3 
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69 
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1 1 
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2 
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691991 

2 

R 

ACT 
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I 67 

7 

18 

1806 
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16 
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16 

90 

691991 

2 

5 

QU 1 

i 3 3 

1 69 

7 

16 

1809 

1 3 

0 

16 

52 

16 

19 

691991 

2 

6 

ROS 

239 

1 69 

7 

18 

18U 

1 3 

0 

18 

15 

19 

2 

691991 

2 

7 

MAD 

169 

1 69 

7 

18 

1019 

1 3 

0 

19 

7 

20 

91 

691991 

2 

8 

4 

MAD 

169 

1 67 

7 

18 

1819 

1 3 

0 

21 

3 

21 

n 

691991 

2 

MAD 

169 

1 69 

7 

18 

1819 

1 3 

0 

21 

10 

21 

39 

691991 

2 

10 

HAD 

166 

1 69 

7 

18 

1 8 1 8 -• 

1 3 

0- 

21 

35 

22 

98 

69 1991" 

2 

11 

WNK 

192 

1 69 

7 

18 

1821 

1 3 

0 

23 

13 

0 

27 

69 |991 

2 

12 

WNK 

192 

1 67 

7 

19 

1821 

1 3 

0 

0 

28 

0 

33 

691991 

2 

1 3 


total no. of files* 30 


PLEASE acknowledge RECEIPj AnD RETURN pORM 22-'Sl - 
DIGITAL’ data ’ ACCOUNT I NQ OpF I C t ’ 

DIGITAL data controller 

NASA GOOOARO SPACE FLIGHT CENTER 

code 565 

GREEnBELT.MD. 20771 ' 

Figure 6.1 , 
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FORMAT OF OGO-F MASTER FILE DATA RECORD 


Word // 

Type 

Variable 


1 

A 

Satellite Name (OGO-F) 

2 

I 

Date (YYMMDD) 

3 

I 

Start Time 

(HHMM) 

4 

I 

Stop Time (HHMM) 

5 

I 

Minutes of 

Coverage 

6 


Not Used 


7 


Not Used 


8 

I 

Pass Number 


9 

I 

Input Tape ; 

Number 

10 


Not Used 


. 11 

I 

Input File ' 

Number 

12 


Not Used 


13 


Not Used 


14 


Not Used 


15 

I 

Date 

Function //I 

16 

1 

User 

Function #1 

17 

A 

Version 

Function //I 

18 

I 

Date 

Function #2 

19 

I 

User 

Function #2 

20 

A 

Version 

Function #2 

21 

I 

Date 

Function //3 

22 

I 

User 

Function //3 

23 

A 

. . Version 

Function #3 

24 

I 

Date 

Function #4 

25 

I 

User 

Function #4 

26 

A 

Version 

Function #4 

27 

I 

Date 

Function //5 

28 

I 

User 

Function //5 

29 

A 

Version 

Function //5 

30 

A 

Film Code 

(CRT/PEN) 

31 

I 

File Status 


32 

I 

Film Number 


33 


Not Used 


34 

I 

10.7 CM Flux Averages 

35 

I 

A Averages 


36 

I 

P wr Sup 

#1 l=0ff 2= On 

37 

I 

Pwr Sup //I ’ 

Time Off /Time On 

38 

I 

HV Pwr Sup 

n l=Off 2=0n 

39 

I 

Pwr Sup iil 

Time Off /Time On 

40 

A 

Scan Code 


41 

I 

Scan Time 


42 


Not Used 



Figure 6.2 
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Word // 


Variable 


Il21 


43 

I 

1st 

Flare 

Number 

44 

I 

1st 

Flare 

Imp Code 

45 

A 

1st 

Flare 

Type Code 

46 

I 

1st 

Flare 

Start Time 

47 

I 

1st 

Flare 

Stop Time 

48 

I 

1st 

Flare 

Max. Time 

49 

I 

2nd 

Flare 

Number 

50 

I 

2nd 

Flare 

Imp Code 

51 

A 

2nd 

Flare 

Type Code 

52 

I 

2nd 

Flare 

Start Time 

53 

I 

2nd 

Flare 

Stop Time 

54 

I ■ 

2nd 

Flare 

Max. Time 

55 

I 

3rd 

Flare 

Number 

56 

I 

3rd 

Flare 

Imp Code 

57 

A 

3rd 

Flare 

Type Code 

58 

I 

3rd 

Flare 

Start Time 

59 

I 

3rd 

Flare 

Stop Time 

60 

I 

3rd 

Flare 

Max. Time 

61 

I 

Start Time Eclipse #1 

62 

I- 

Duration 

Eclipse #1 

63 

I 

Start Time Eclipse #2 

64 

■ I 

Duration 

Eclipse #2 

65 

I 

Start Time Eclipse //3 

66 

I 

Duration 

Eclipse #3 


. WORDS 67-100 ARE NOT USED 


Figure 6.2 
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FORMAT OF PUNCHED CARD GENERATED FROM MASTER FILE DATA RECORD 


Columns Code Description 


1-4 

A4 

Satellite Name 

6-11 

16 

Date (YYMMDD) 

13-16 

14 

Start Time (HHMM) 

18-21 

14 

Stop Time (HHMM) 

26-30 

15 

Input Tape Number 

32-33 

12 

Input File Number 

76-77 

12 

User Code 

79-80 

12 

Function Code (1) 


Figure 6.3 
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Figure 6.4 
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FORMAT & DESCRIPTION OF MASTER FILE PRINTOUT 


Print Positions 


Description 


2-7 

10-13 

17-20 

23-25 

28-31 

34-35 

39-42 

45-46 

49 

52-55 

58 

61-64 

67-69 

72-75 

78-82 

84 

85 

88-91 

93-96 

98-101 

105-108 

112-113 

118-121 

125-126 


Date of Acquisition (YYMMDD) 
Start Time of Acquisition (HHMM) 
Stop Time of Acquisition (HHMM) 
Minutes of Coverage (MMM) 

Input Tape Number 

Input File Number 

Daily 10,7 CM Flux Average 

Daily Ap Averages 

HV PWR Sup #1 l=0ff, 2=0n 

PPJR Sup #1 Time Off/Time On 

HV PWR Sup //2 l=0ff, 2=0n 

PWR Sup #2 Time Off /Time On 

Scan Code 

Scan Time 

1st Flare Number 

1st Flare Imp Code 

1st Flare Type Code 

1st Flare Start Time 

1st Flare Stop Time 

1st Flare Max, Time 

Start Time Eclipse #1 

Duration (MM) Eclipse //I 

Start Time Eclipse //2 

Duration (MM) Eclipse #2 


Figure 6 .5 
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PUNCHED 

CARD FORMAT 

TO UPDATE PRODUCTION MASTER FILE 

Columns 

Code 

Description 

1-4 

A4 

Satellite Name 

6-11 

16 

Date (YYMMDD) 

13-16 

14 

Start Time (HHMM) 

18-21 

14 

Stop Time (HHMM) 

26-30 

15 

Tape Number Processed 

32-33 

12 

File Number Processed 

38 

11 

Plot Code jl)=No Plot 
J 1=CRT Plot 
^=PEN Plot 

40-41 

12 

File Status (see Figure 6.7) 

45-46 

12 

Film Number (CRT only) 

76-77 

12 

User Code 

79-80 

12 

Function Code (1) 


Figure 6,6 
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FILE STATUS - NUMERICAL CODE DEFINITIONS 


0 

1 

2 

3 

4 

5 

6 

7 

8 


ACCEPTABLE FILE 

INSUFFICIENT DATA IN FILE (HEADER ONLY) 
DAMAGED DATA TAPE 
POWER SUPPLIES OFF 

FILE CONTAINS LONG AND SHORT SCAN DATA 
MISSING DATA TAPE 

DISCREPANCY BETWEEN START AND STOP TIMES 
BAD DATA IN FILE 

FOUR (4) MESSAGES PER SCAN POSITION 


Figure 6.7 
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7.0 OGO-F INTERACTIVE GRAPHICS APPLICATION PROGRAM (OGQGRAF) 


Program OGOGRAF is an interactive graphics program designed to 
aid in the analysis of spectral data acquired from the six sensors 
flown on the OGO-F satellite. The batch processing program^ RDLP, 
which performed bulk analysis of the digitized data was modified to 
create a data base for OGOGRAF, the graphics application program. 

A format was developed to include all possible data for graphics 
analysis. The OGO-F data on the original digitized data tapes was 
comprised of multiple files, one file for each acquisition. Each 
acquisition consisted of several scans of data. Each scan contained 
512 samples for each of six sensors. Each sensor recorded data from 
a different spectral range. 

A file of data for the graphics application data base was 
structured as follows. A file header record contained the original 
tape and file number, date of acquisition, orbit number and 20 
parameters used in the batch processing program. The 20 parameters 
were included since they were control card inputs and could be 
changed in the batch processing program. These parameters were used 
primarily in the aspect and backlash correction routines. Each scan 
in the file was preceded by a scan header record which included the 
niamber of data records in the scan and whether the scan was a back- 
ward or forward scan. Each data record contained the scan position, 
time, wavelength, uncorrected sensor count, and aspect corrected 
sensor count for each of the six sensors and also the two aspect 
angles with the three components of each angle. This format and file 
structure included all parameters necessary for any stage of analysis 
from raw data to corrected data. 

The graphics application program was designed to anticipate a 
multiple file data base. The first interactive task was for the 
researcher to enter via the alphanumeric keyboard the tape and file 
number of the data set to be analyzed. The character string was 
decoded and edited. If an illegal character was decoded, an explana- 
tory message was displayed. If the entire data base was searched 
without encountering the requested file, a prompter message was 
displayed on the screen requesting a new file or program termination. 
Once the file was located, the researcher was requested to define, 
via the keyboard, the log scaling desired for the initial display. 

This scaling could be changed any time during the analysis by picking 
the Y-scale button and entering the new Y-maximum and minimum through 
the keyboard. 

The data for the first scan of the file was then stored in the 
random access file. All data in each record were stored to facilitate 
the displaying of other parameters from this scan at the researcher’s 
discretion. As the data were stored, the wavelengths and corrected 
sensor values for the first of the six sensors were stored in the 
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display arrays. As they were being stored, the minimum and maximum 
wavelengths were saved. These wavelengths were used to define X- 
minimum and X-maximum for the console display. This allowed for the 
most efficient utilization of the limited expanse of the visible 
console screen. The sensor values were converted to logs, and the 
logs and the wavelengths were converted to graphics display units. 

The display itself consisted of the axes, labeling of the axes, start 
and stop time of the scan, the scan direction and the 512 samples 
connected by line segments. Figure 7.1 is a plot of this initial 
display. 

The primary functions of OGOGRAF were five in number. These 
functions were display flexibility, peak analysis, experiment back- 
ground definition and application, sensor correction analysis and 
documentation of interactive graphics procedures and results. 

Display flexibility included the ability to advance from scan to 
scan, to select the sensor within a scan being displayed, to choose 
the parameters for a particular sensor to be displayed and to modify 
both the X and Y scaling of the display. 

The NEWSCAN primary button, when picked, automatically read the 
next consecutive scan, stored the data in the random access file and 
displayed the first of the six sensors. The SENSOR primary button, 
when picked, displayed a menu of the six sensors with appropriate 
pseudonyms. The researcher then picked the sensor to be displayed. 

OGOGRAF then read the random access file extracting the sensor readings 
and wavelengths for that sensor, converted the sensor readings to 
logs, converted the logs and wavelengths to display units and formed 
the display with appropriate labels. Thus the researcher could dis- 
play the data from any of the six sensors in any order. If he 
decided to display the sensors in a predefined sequence, the picking 
of the CONTINUE primary button advanced to the next of the six 
sensors and displayed the data. The sensor sequence was 3600, 2160, 

1800, 1500, 1200 and 900 grating lines/MM. Normally, the aspect 
corrected sensor values were displayed. By picking the RAW DATA 
primary button, the uncorrected sensor data would be extracted from 
the random access file. Any future requests for new sensors would 
display the uncorrected data for that sensor. Figure 7.2 is a plot 
of the uncorrected data. The COMPUTED primary button resulted in the 
corrected sensor data being extracted and displayed. The Y-SCALE primary 
button allowed the currently displayed data to be rescaled on the 
log axis. The minimum and maximum decade values were entered via the 
alphanumeric keyboard. The newly defined Y-Scaling would also apply 
to any new displays until Y-SCALE was picked again, and new log scaling 
was defined. Figure 7.3 is a plot of the data in Figure 7.1 with modified 
scaling. The SCISS0RS primary button was used to select any segment 
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of the current display and to redisplay that segment over the full 
extent of the X-AXIS. The researcher defined the leftmost and right- 
most extent of the segment by moving the tracking cross to the desired 
X locations. Two vertical lines were displayed at the extents. If the 
researcher desired that segment, he picked the ACCEPT secondary button; 
and the data in the display arrays within the boundaries of the segment 
were extracted, a new X-minimum and maximum defined; and the data were 
displayed over the entire X-AXIS of the display. This allowed the 
researcher to more clearly analyze portions of the spectrum with a 
high density of lines. The SCISS0RS application applied only to the 
data currently being displayed. A request for a new sensor or a new 
scan would scale the entire data set for the sensor along the X-AXIS. 

Figure 7.4 is a plot of the first half of Figure 7.1 expanded. 

By selecting the appropriate task request buttons which pertain 
to display manipulation, the researcher had full flexibility to display 
any portion of any sensor of any scan in his data set. A new file in 
the graphics database could be accessed by selecting the PINIT primary 
button which accepted the new tape and file number, located it in the 
data base, stored the first scan in the random access file and displayed 
the data from the 3600 sensor on the first scan. 

The second primary function of OGOGRAF was to allow the researcher to 
interactively define peaks from the data currently being displayed. The 
VALUES primary button allowed the researcher to position the tracking 
cross over a data sample. The ACCEPT secondary button was picked, and the 
data sample closest to the center of the tracking cross was flashed. The 
sensor value and the wavelength of the picked point were also displayed 
on the screen. By selecting the NEXT secondary button the researcher 
could display the coordinates of as many data samples as desired. The 
task was completed when the END secondary button was chosen. The PEAK 
primary button allowed the researcher to pick data samples considered to 
best represent peaks for the sensor being displayed. Since most scans 
contained 20-30 legitimate, identifiable peaks, it was necessary to 
simplify this procedure as much as possible. Consequently, when PEAK 
was chosen, the sensor readings being displayed were interrogated. Those 
which fulfilled a simple possible peak criteria were made light -pen 
sensitive on the display. The data samples not considered to. be a possible 
peak remained insensitive to the light-pen. Forty-fifty samples of the 
512 fulfilled the peak criteria. Most importantly, the light-pen sensitive 
samples would be bordered by insensitive samples to greatly reduce the 
possibility of picking an unwanted sample as a peak. By displaying the 
possible peaks as light -pen sensitive points, it was only necessary for 
the researcher to engage the switch on the light-pen to return the iden- 
tity of the sample to the processing module in the CDC 6600. The 
researcher knew the sample picked because it flashed. If it was hot 
the correct sample, it would be deleted from the queue by simply again 
picking the flashing point on the display. In this manner, he could 
pick all peaks in a display in a matter of seconds, review the flashing 
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points, make any additions or deletions, and pick the END secondary 
button. When END was picked, the queue of selected peaks was trans- 
ferred to the 6600. The time, wavelength and intensity of the peaks 
were extracted from the random access file and printed for a 
permanent, hard copy record of the analysis. 

The third primary function of OGOGRAF was to analyze and to define 
background intensities to study the consistency of the background and the 
subtraction of the background from peak values to obtain the true peak 
value. Each of the six sensors had a unique background profile. Some 
had a constant background throughout the scan, others had segments with 
constant background and segments with exponentially increasing background. 

Two primary buttons were made available for this analysis. BKGROUND 
enabled the researcher to divide the scan being displayed with the 
tracking cross into one, two or three segments which accommodated all 
background combinations for the six sensors. Vertical lines were dis- 
played at the segment interfaces. If the lines showed that the segments 
were not correct, they could be redefined via the tracking cross; and 
new vertical lines would be displayed. The type of background, constant 
or exponential, and the actual backgrounds were then defined for each 
segment. The researcher picked either the HORIZONTAL or EXPONENTIAL 
secondary buttons with the light-pen. For the former, he moved the 
tracking cross to the desired Y-Value to be used as the background for 
that segment. He picked ACCEPT, and a horizontal line was displayed 
over the range of that segment. If the background looked too high or 
too low, it could be redefined with the tracking cross. When ACCEPT 
was picked a second time, the segment boundaries in angstroms and 
the background level were stored in an array for that sensor. For 
the exponential background, the tracking cross was used to define two 
points in the segment. Since the display was semi-log, an exponential 
background line could be defined by a straight line on the log scale. 

The line defined by the two points was displayed over the extent of 

the segment. If acceptable, the segment boundaries in angstroms and the 

slope and Y- Intercept of the line were saved in the array for that sensor. 

If the background line was unacceptable, two new points could be defined 
with the tracking cross. Figure 7.5 shows the 1200 lines/MM sensor before 
background subtraction. 

Once the background lines had been defined for a sensor, the re- 
searcher could subtract this background from all 512 samples in the scan 
by picking the SUBTRACT primary button. The wavelength for each sample 
would be matched with the applicable background segment. If it was a 
constant background segment, the background defined by the horizontal 
line would be subtracted. If it was an exponential background segment, 
the first degree equation would be solved for the wavelength yielding 
the log of the background. The antilogs of the sensor reading and 
background were calculated; the background was subtracted; and the reading 
minus the background would be converted to a Log^^ value. When all 512 samples 
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were processed, the new scan minus the background was redisplayed. By 
rescaling the X-AXIS, the researcher could then analyze the results. 


Once the background parameters had been established for a sensor, the 
data for that sensor for any other scan could be corrected for background 
by selecting the SUBTRACT primary button after the new sensor data had 
been displayed. Analysis of the consistency of the backgound levels was 
performed in this manner. However, the background parameters could be 
redefined at any time by selecting the BKGROUND primary button. 

The fourth primary function of OGOGRAF was the sensor correction analy- 
ses. During the bulk processing of the experiment data, corrections 
were made for angular aspect, sensor efficiency and mechanical backlash 
when the 512 step collimator changed direction. Also, the determination 
of the spectral range was inconsistent between forward and backward 
scans. Twenty parameters used in these correction and wavelength de- 
termination formulas were control card inputs. The CNST primary button 
when selected displayed the values of these 20 constants. The researcher 
could change one or more of them via the alphanumeric keyboard. If the 
COMPUTED primary button were then selected, it displayed another option 
secondary button, RE-COMPUTE. If this button were selected, the raw 
data for the sensor being displayed would be extracted from the random 
access file. Sensor correction for angular aspect would be recalculated 
along with the wavelengths using the same processing modules as in the 
batch processing program except with the newly defined constant (s) . The 
new corrected sensor values were then displayed against the new wavelength 
range. Figure 7.6 shows the data from Figure 7.1 recalculated with 
a redefined backlash correction. Note that the wavelength scaling has 
been modified. 

The fifth primary function of OGOGRAF was the ability to document 
procedures and results. In many instances the researcher could not 
compare results of analysis between scans or between displays before and 
after subtle changes in the correction parameters. The PLOT button , 
allowed the researcher to create a 35MM film frame of the display at / 
that moment. A 50-character message could be entered via the keyboard. 

This allowed for a permanent record of all pertinent stages of analysis 
and processing which could be further analyzed in the future or compared 
to results from prior graphics sessions. The figures used in this 
section are plots generated by the PLOT primary button. 

The execution of program OGOGRAF requires that the file or files to 
be processed be stored as permanent disk files. Graphics application 
programs cannot request magnetic tape input. Following is the deck setup 
for processing a permanent file cataloged as GRAFDATA. Note that 
Graphics Programs must be submitted to be run on system A. 
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Card 


1 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 



Orange 6600 Input Card (System A) 

Core Memory = 77K 
C.P. = 200 seconds 

P.P. = 500 seconds 

Plot 

Write "Graphics Job" in special instructions box 
DELOl, CM77777, T200. 7276 DELANEY 

ATTACH, TAPEl, GRAFDATA, MR=1. 

ATTACH, CRT, NEWONLINECRT , MR=1. 

FTN, A, 0PT=2. 

CDPYBR, CRT, LGO, 27. 

FTN, A, 0PT=2. 

LOADER, PPLOADR. 

NOREDUCE. 

LOAD, LGO. 

NOGO 

NOREDUCE. 

AEFILE, JTl, SOURCE. 

SOURCE. 

RELEASE, JTl. 

EXIT. 

EOR (7/8/9 in column 1) 

THE OGOGRAF FORTRAN CARD DECK. 


EOR 

JTl 

EOR 

EOF (6/7/8Z9) 

N.B. All cards with the exception of those in the Fortran card 
deck begin in column 1. Punctuation must be followed. 
Interactive processing can begin when the tape and file 
prompter messages appear on the graphics console screen. 
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